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Abstract 

Based on a successful description of II- VI ternary alloys, which introduces an 
empirical bowing parameter to the widely used virtual crystal approximation, 
we set up a tight-binding Hamiltonian to describe the Zni_j^CdySei_a;Tea; 
and Zn.9Cd.1S.07Se. 93 quaternary alloys. We just use a formula that can be 
thought as a straightforward generalization of the virtual crystal approxima- 
tion for this case. Our Hamiltonians reproduce very well the change in the 

band gap value with the composition observed in recent experimental reports. 
PACS: 71.15.Fv, 71.55.Gs, 73.10.-r 
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Typeset using REVT^ 



Interest in optical devices that can operate in the visible spectrum range has motivated 
a vigorous research in II-VI wide band gap semiconductor compounds. The main goal in 
these studies is to seek for the independent control of the fundamental parameters of the 
compound: the band gap, the lattice constant and the valence band offset. The interest 
has been recently focused on the study of the ternary and quaternary alloys. In the quater- 
nary alloys it is the simultaneous substitution of the anion and the cation that allows the 
manufacture of a compound with band gap and the lattice constant values within certain 
limits. This procedure is very important in optoelectronic engineering. In fact, it could 
be produced a binary/quaternary heteroj unction with no stress at the interface and with a 
predetermined value for the band offset fTj-Q. 

Recently, Brasil et al. have reported photoconductivity and photoluminiscence exper- 
iments on the Zni_j^CdySei_2:Tej; system, and they have fitted the measured band gap value 
to a quadratic function of the composition; Ichino et al. have used the Zn.9Cd.1S.07Se. 93 
quaternary alloy as an active layer in a laser diode operating in the blue-green to the UV 
spectral region. 

In a previous work, we have described the change of the band gap value for the ternary 
alloys and compared the results with the existing experimental data To describe the 
ternary alloys we have used carefully studied tight-binding Hamiltonians for the binary 
compounds 0. We took into account nearest neighbors interactions with an orthogonal 
basis of five atomic orbitals {sp^s*) per atom; in this approximation the s*-state is included 
to better describe the lower states in the conduction band. We have used the two center 
integrals approximation in the Slater-Koster language and we took into account the spin- 
orbit interaction 0,0]. The tight-binding method describes rather well the valence band 
region and gives correctly the band gap value ^j. To include composition we have generalized 
the virtual crystal approximation (VGA) by introducing an empirical bowing parameter. 
As we have found, the VGA model can describe properly the behavior of the band gap 
if we introduce an empirical bowing parameter to calculate the "s" on-site tight-binding 
parameter of the substituted ion @]. 
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In the semiconductor compounds the band gap, Eg, is given by the difference between 
the minimum of the conduction band, E^, and the maximum of the valence band, E^, 

Eg = Ec — Ey (1) 

In particular, for the compounds with direct transition the expressions for E^. and i^^^ at F , 
in the tight-binding method, are: 



E. = + Ji^^^y + Vi (2) 



and 



r-, Epa + Epc + Aa + Ac / / Epa — Epc + Xa — Ac\2 

Ev = Y ) + Vxx l-^j 

where Eav = s, p; v = a, c) are the on-site tight-binding parameters with s or p 
character for the anion (a) and the cation (c), Vaa are the nearest neighbor interaction 
parameters, and \y the spin-orbit parameters. 

To describe the alloy we have used the VGA model, as we stated above. In this approxi- 
mation the tight-binding parameters are given by the weighted averages of the corresponding 
end-point parameters. For a ternary alloy ABi_a;Cx-, for example, we have: 

E^c.'{x)=xEt2> + {l-x)Et^, (4) 

where E^ is the corresponding tight-binding parameter of the different binary compounds 
and a, a' denote the atomic orbitals {sp^s*). We will use this expression for all the tight- 
binding parameter but the "s" on-site ones. 

As it is well known, the VGA model by itself does not describe correctly the non-linear 
behavior of the band gap value in the alloy. Therefore, to take into account the non-linear 
behavior of the band gap in the ternary description, we have included an empirical bowing 
parameter in the "s" on-site tight-binding parameter of the substituted ion. This is because 
the "s" on-site tight-binding parameter is responsible for the correct energy position of the 
conduction band in the F-point IQ, as can be inferred from eq. (2). In an explicit way. 



when we have an anion substitution we introduce an empirical bowing parameter in the "s" 
on-site anion- VGA expression, and the same for the cation case. We use the next VGA 
expression for the "s" on-site tight-binding parameter 

E'^^{x,bu) = Es^{x) + bux{l - x), u = a, c (5) 

where Esv{x) is given by eq. and is the empirical bowing parameter appropriate to 
each substitution. We define the empirical bowing parameter as 



ss ss 



where El'if^ is the "s" on-site tight-binding parameter for the compound 1 (2), V^^*-^-* the 
corresponding nearest neighbour s — s interaction parameter for the compound 1 (2), and 
we use the sign +(— ) for cation (anion) substitution, the proportionality constant is taken 
as k = 1/8, and we take A as our free parameter in order to obtain the best fit to the 
experimental data. For the most of the alloys we have used A = 1.75 0], and for the S- 
based alloys we have used A = 1.0. We will show that these values for A give us good 
agreement with experiment. The values for the empirical bowing parameter, given in Table 
I, are the same that we have used in the ternary description (see Ref. [Q). With this 
Hamiltonians we can calculate the band gap for any composition in a straightforward way. 
However, before to do that we should to extend our approximation in order to incorporate 
properly the quaternary alloys. 

In the quaternary alloy Ai^yByCi-x^x the VGA expression for the tight-binding param- 
eters, that is the appropriate extension of eq. (^), is 

E^a'ix, y) = xyEl^, + x(l - y)E^^, + (1 - + (1 - a;)(l - y)E^2'^ (7) 

where the E'-^'s are defined in eq. (^). In these alloys we have both substitutions, anion- as 
well as cation-substitution. For that reason we have found natural to extend (in the spirit 
of the VGA) our formula, for the "s" on-site parameters, of the ternary alloys as follows: 
For the anion substitution, 



Esa{x, y, ha) = Esa{x, y) + x(l - + (1 - y)h^), (8) 

where Esa{x, y) is given by eq. (0), 6^ and are the empirical bowing parameters for the 
anion substitution in the BCi_a;Da; and ACi-ajD^; systems, respectively. 
For the cation substitution, 

y, K) = Escix, y) + y{l - + (1 - (9) 

where Esc{x,y) is given by eq. (|^), and 6^ are the empirical bowing parameters for 
the cation substitution in the Ai_yByD and Ai_yByC systems, respectively. Note that in 
the end- values of the compositional variable y in eq. (or x in eq. (^) we obtain the 
appropriate expression for the ternary case, eq. (5). Using these expressions, for the "s" 
on-site tight-binding parameters, we can calculate the band gap value and the electronic 
structure for any composition of the quaternary alloys. 

Table II shows a comparison of our calculation of Eg{x,y) and the photoconductivity 
measurements of Brasil et al. for Zni-yCdySei-^^Tei, As can be judged through the study 
of the table our calculation agrees well with the measured values. In general, we reproduce 
the experimental values within 2% accuracy. 

Fig. 1 shows a 3-D graph for Eg{x,y). The dots are the experimental data, the mesh 
is our calculation. Note that all the experimental data are just over the calculate surface. 
The edges are the four ternary boundaries. As we have noticed, our calculation reproduces 
well the known bowing of Eg, both in the anion substitution as in the cation case. In the 
anion substitution the bowing of Eg is more noticeably than in the cation one. This fact 
can be inferred from the figure if we take the projection of Eg on the X — Z plane. In 
the same way, we also appreciate the cuasi-linear behavior of Eg in the cation substitution 
case. We can see this facts from the Table I as well. Notice that, in absolute value, the 
bowing parameter for the anion substitution is greater than the cation one. Although our 
empirical bowing parameter is not given by the experiment it can be shown, through the 
tight-binding equations, that a great empirical bowing parameter, in absolute value, gives 



5 



great optical bowing parameter; and the optical bowing parameter is the measured value in 
the experiment (see Ref. 0] for details in this fact). 

We have calculated the electronic structure for the Zn.gCd.i 8.0786.93 quaternary alloy, as 
well. The particular interest for this quaternary alloy is the experimental report of Ichino 
et al. These authors propose the alloy as an active layer in laser diode operating in the 
blue-green to ultraviolet spectral region. Fig. 2 shows our calculated electronic structure. 
This electronic structure is representative of a semiconductor compound with direct band 
gap. From the figure we appreciate that the lowest conduction band shows lesser dispersion 
than the calculated one, using the empirical pseudopotential method (EPM), by Feng et al. 
0. However, the general pattern obtained for the valence band is the expected one from 
tight-binding calculations. The calculated band gap for this alloy is 2.705 eV that is in good 
agreement with the experimental value of 2.73 eV at 4.2 K given by Ichino et al. 0, and 
better than the value of 2.648 eV calculated by Feng et al. using the EPM. 

It is noteworthy that: This work is based on a correct description for the ternary systems. 
We believe that this fact gives an important support to our method in general and to the 
potential use of this Hamiltonians in further calculations for surfaces, interfaces, quantum 
wells and superlattices of II- VI compound systems . 

In conclusion, we have calculated the changes in the band gap value of the 
Zni_j^Cdy8ei_^.Te2: and Zu gCd.i 8.0786.93 quaternary alloys within a tight-binding descrip- 
tion. For the pure binary compounds, we have used the tight-binding parameters which 
describe well the known band structures. To describe the alloys, we have used the virtual 
crystal approximation reformulated according to reproduce the observed non linear behavior 
of the band gap with the composition. We have introduced an empirical bowing parameter 
in the tight-binding parameters of the "s" atomic orbitals which are known to be responsible 
for the correct energy position of the conduction band in the F-point. Our results agree well 
with experimental results. Our Hamiltonians can be used as a basis for other calculations 
that include this kind of pseudobinary compounds. 
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TABLES 

TABLE I. Empirical bowing parameters obtained using the eq. 
parameters given in the Refs. [5,11] 



(6) and the tight-binding 



Compound 




b„ 


be 


ZnSei_a;Te2; 




-6.964 


- 


CdSei_a;Tea; 




-0.195 




CdSi-o^Se^ 




-0.136 




ZnSi_a;Se3; 




-2.833 




Zni_a;Cda;Se 






0.037 


Zni_^Cd^Te 






0.020 


Zni_j;Cdj;Se 






1.349 


TABLE IL Measured values for the band 


gap given in Ref. [1] for the Zni_ 


yCdySei_j;Tej; 


quaternary alloys 


are compared with our calculation. Band gap values given in eV. 




X 


y 




£;g(theo)'^ 


0.015 


0.10 


2.610 


2.642 


0.020 


0.23 




2.453 


0.050 


0.22 


2.443 


2.427 


0.070 


0.34 


2.244 


2.256 


0.080 


0.28 


2.362 


2.316 


0.038 


0.32 


2.024 


2.002 


0.039 


0.11 


2.224 


2.194 


1.000 


0.31 


2.049 


2.099 



^ Values taken from photoconductivity measurements of Brasil et al. Ref. [1] 



^ This work 
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FIGURES 

FIG. 1. 3-D representation for the calculated band gap of the Zni_j^Cdj,Sei_a;Tea; quaternary 

alloy, as function of the composition {x,y), using the tight binding method and the virtual crystal 
approximation as is proposed in this work. The points are the experimental data of Brasil et al. 
[1]- 

FIG. 2. Electronic structure for the Zn.9Cd.1S.07Se.93 quaternary alloy calculated using our 
tight-binding parametrization. 
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